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The origin of unconventional superconductivity, including high-temperature 
and heavy-fermion superconductivity, is still a matter of controversy. Spin ex- 
citations instead of phonons are thought to be responsible for the formation of 
Cooper pairs. Using inelastic neutron scattering, we present the first in-depth 
study of the magnetic excitation spectrum in momentum and energy space in 
the superconducting and the normal states of CeCu2Si2- A clear spin excita- 
tion gap is observed in the superconducting state. We determine a lowering of 
the magnetic exchange energy in the superconducting state, in an amount con- 
siderably larger than the superconducting condensation energy. Our findings 
identify the antiferromagnetic excitations as the major driving force for super- 
conducting pairing in this prototypical heavy-fermion compound located near 
an antiferromagnetic quantum critical point. 

While conventional superconductivity (SC) is generally incompatible with magnetism, 
magnetic excitations seem to play an important role in the Cooper pair formation of un- 
conventional superconductors such as the high-T c cuprates or the low-T c organic and heavy- 
fermion (HF) superconductors. Since the discovery of SC in CeC^Si^, antiferromagnetic 
(AF) spin excitations have been proposed as a viable mechanism for SC£~-. The discovery 
of SC at the boundary of AF order in CePd2Si^ has pushed this notion into the framework 
of AF quantum criticality-. Unfortunately, such quantum critical points (QCPs) proximate 
to HF superconductors typically arise under pressure, which makes it difficult to probe their 
magnetic excitation spectrum. 

Here, we report a detailed study of the magnetic excitations in CeCu2Si2, which exhibits 
SC below T c pa 0.6 K. This prototypical HF compound is ideally suited for our purpose, 
since SC here is in proximity to an AF QCP already at ambient pressure (cf. Fig. Ufa)). 
As displayed in Fig. [T^b) CeCu 2 Si 2 crystallises in a structure with body-centred tetrago- 
nal symmetry and is one of the best studied HF superconductors and well characterised 
by low-temperature transport and thermodynamic measurements^. Moreover, those mea- 
surements in the field-induced normal state have already provided evidence that the QCP 
in this compound is of the three-dimensional (3D) spin-density-wave (SDW) type^. The 
spatial anisotropy of the spin fluctuations in superconducting CeCu 2 Si 2 was measured at 
T = 0.06 K and at an energy transfer Tvjj = 0.2 meV and is shown in Fig. QJc). These mag- 
netic correlations display only a small anisotropy (a factor of 1.5) in the correlation lengths 



2 



between the [110] and the [001] direction. Therefore, these quite isotropic spin fluctuations 
are in line with thermodynamic and transport measurements exhibiting C /T = 70 — ay/T 
or p — po = AT a , a = 1 — 1.5&2, and strongly support a three-dimensional quantum critical 
SDW scenario^ 4 . We are able to identify the magnetic excitations in the normal state of 
paramagnetic, superconducting CeCu2Si2, around the incommensurate wave vector- of the 
SDW order nearby in the phase diagram (cf. Fig. QJa)), and further establish the system's 
proximity to the AF QCP through the observation of a considerable slowing down in the 
spin dynamics. Going into the superconducting state, a spin gap opens out of a broadened 
quasielastic response which extends to much higher frequencies (10 x the superconducting 
gap). These data allow us to establish a saving in the AF exchange energy that is con- 
siderably larger than the superconducting condensation energy, thereby providing the first 
demonstration of the nearly-quantum-critical AF excitations as the major driving force for 
unconventional SC. 

Superconductivity and antiferromagnetism in CeCu2Si2 

The SC in CeCu2Si2 we consider is close to the AF QCP, and is to be contrasted with a 
second superconducting dome appearing at high pressure which is thought to be associated 
with a valence instability and the concomitant fluctuations*^. This AF QCP is located 
within the narrow homogeneity range of the "122" phase in the ternary chemical Ce-Cu-Si 
phase diagram of this tetragonal compound^. Correspondingly, we can prepare homogeneous 
samples (with slight Cu deficit) from the antiferromagnetically ordered side (A-type) and 
(with tiny Cu excess) from the paramagnetic, superconducting side of the QCP (S-type); by 
contrast, crystals very close to the 1:2:2 stoichiometry exhibit a ground state where SC and 
AF compete with each other without microscopic coexistence (A/S-type)^. The AF order 
was found to be an incommensurate SDW&. At T = 0.05 K, well below T N « 0.8 K, the A- 
type CeCu2Si2 exhibits an ordered magnetic moment p or( i « 0.1 p-Q and an incommensurate 
propagation vector r ps (0.215 0.215 0.53). The latter can be ascribed to a nesting wave 
vector of the renormalised Fermi surface. However, a full microscopic description of the 
magnetic order remains to be addressed. 

In order to study the superconducting state in CeCu2Si2 in detail, neutron scattering 
results^ on an S-type single crystal are presented in this Article. A previous experiment 
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was severely hampered by a large experimental background and a low signal-to-background 
ratio^. Thermodynamic measurements confirmed that this crystal is superconducting with 
a T c w 0.6 K and an upper critical field B c2 {T = 0) < 2T— . Elastic neutron scattering 
measurements did not feature resolution-limited magnetic Bragg peaks in S-type CeCu2Si2 
in accordance with thermodynamic measurements. However, at positions where magnetic 
satellite peaks are observed in A-type CeCh^Si^, e.g., at Qaf = (0.215 0.215 1.458), relative 
to a nearby nuclear Bragg reflection G (Qaf = G±r), the S-type crystal exhibits quite weak 
correlation peaks at low temperatures^. They are still present above T c and disappear at 
T m 0.8 K, very similar to the behaviour of the SDW order in A-type CeC^Si^. Although 
these peaks were found to be purely elastic within the energy resolution (« 57/xeV), their 
linewidth in Q space is considerably broadened corresponding to a correlation length of 
50 — 60 A being comparable to the superconducting coherence length of order 100 A-. Thus, 
static magnetically ordered regions seem to exist in a quite small part of the sample and are 
separated from the surrounding superconducting regions. 

Spin dynamics in CeCu2Si2 

We probe the magnetic response of CeCu2Si2 through extensive inelastic neutron scat- 
tering measurements around Q = Qaf = (0.215 0.215 1.458), since no appreciable magnetic 
intensity has been detected elsewhere in the Brillouin zone. Fig. 2(a) displays energy scans 
at this Qaf position and at a general position Q = Q arb = (0.1 0.1 1.6), where no cor- 
relation peaks emerge, but which has the same |Q|. Both data sets were recorded in the 
superconducting state at T = 0.07 K. At Q ar b only the incoherent elastic background con- 
tribution with instrument resolution is seen, while no magnetic intensity could be detected. 
In contrast, at Qaf the response shows a strong inelastic signal with a long tail of the in- 
tensity extending beyond fiw = 2meV (cf. inset of Fig. 2(a)). The missing spectral weight 
at low energies is an indication for a spin-excitation gap in the superconducting state. The 
spectrum recovers the missing weight at the gap edge, thereby constituting an inelastic 
line. The data can be described by a quasielastic Lorentzian line with a spin excitation gap 
fi^gap ~ 0.2 me V and with a density of states as for the electronic gap of a d-wave BCS 
superconductor (solid lines in Figs. 2(a) and 3(a)). ftcjgap ~ 3.9k B T c is found to be 10% 
smaller than the value predicted for a weak-coupling d-wave superconductor— and falls 20% 
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below 2A //c_bT c = 5.0 as determined by Cu-NQR for CeCi^Si^ 1 ^. To unambiguously 
relate the inelastic magnetic excitation to the superconducting state, it was necessary to 
perform additional measurements in the normal state. 

Energy scans recorded at Qaf in the normal state are shown in Fig. 2(b). Notably, 
independent of how the normal state is reached, i.e., above T c at T = 0.8 K and B = 
or above B& at T = 0.07 K and B = 2T, the magnetic response is almost identical and 
appears to be quasielastic. The fits to the quasielastic magnetic response with a Lorentzian 
lineshape give a good description of the data as seen in Figs. 2(b) and 3(a). With increasing 
temperature in the normal state the magnetic response weakens in intensity and broadens 
considerably. Starting from T ps 0.11 meV at T = 0.07 K the linewidth of the quasielastic 
response at Qaf increases to T w 0.235 meV at T = 1.7K (Fig. 3(b)). This considerable 
slowing down of the response when lowering the temperature indicates the proximity of 
S-type CeCu2Si2 to the AF QCP. T(T) extrapolates to a finite value at T — > 0, since the 
S-type single crystal is located on the paramagnetic side of the QCP (cf. Fig. 1(a)). A 
related critical slowing down was observed in magnetically ordered A- type CeC^Si^. 

The fact that the magnetic excitation gap disappears in the normal state, i.e., above T c , 
and also above B c2 at low temperatures, where the magnetic short-range correlations still 
persist, gives direct evidence that the spin gap huj ga , p is related to the superconducting state. 
Its temperature variation is displayed in Fig. 3(c) and has been derived from fits to the data 
shown in Figs. 2(a) and 3(a) and additional scans. As indicated by the solid line, hu) gap 
follows, within the error bars, the BCS form for the superconducting gap amplitude 2A(T). 

We now turn to the momentum dependence of the magnetic response around Qaf in the 
superconducting state. Fig. 4 displays Q scans along (h h 1.458) across Qaf recorded at 
different energy transfers hu and at T = 0.06 K. The single peak seen at hu w 0.2 meV splits 
upon increasing energy transfer into two peaks which move further apart from each other 
along with a marked decrease in intensity. Fits with two peaks of Gaussian lineshape (solid 
lines) yield a good description of the data. The peak positions for different ftw, drawn in the 
inset of Fig. 4(a), yield a linear dispersion relation. We conclude that the spin excitations 
are part of an overdamped dispersive mode. Its velocity as read off the slope of the dispersion 
curve, v exc = (AAA ± 0.86) meVA, is substantially smaller than the strongly renormalised 
Fermi velocity ~ 57 me V A- (lmeVA = 153 m/s). This indicates a retardation of the 
coupling between the heavy quasiparticles and the quantum-critical spin excitations. 
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Superconducting condensation and magnetic exchange energies 

The observed spin excitations both below and above T c allow us to estimate the decrease 
of magnetic exchange energy in the superconducting state as compared to the putative 
normal state. This saving of exchange energy is determined as follows^i^: 



where E x is respectively the exchange energy in the normal (N) and superconducting (S) 
states, <> indicates an average over the first Brillouin zone, and q = (q x ,Qy,Qz) denotes 
a momentum transfer in the first Brillouin zone, i.e., Q = G + q. /(q) is the exchange 
interaction between the localised /-moments and contains nearest (ii) and next nearest (I2) 
neighbour terms: 



where a and c are the lattice constants, and the precise form of /2 is given in^ 4 -. The 
inclusion of the next nearest neighbour terms is a consequence of the three-dimensional 
nature of the spin excitations of CeCu2Si2 (cf. Fig. CE^c)). This is different from the cuprate 
superconductors and e.g. CeCoIns, where the observed behaviour is predominantly two- 
dimensional. As described in detail in the supplementary material^ 4 -, we find a magnetic 
exchange energy saving of AE X = r\ 4.8 • 10 _3 meV per Ce (77 ~ 1.25, 77 being a measure of the 
SC volume fraction^ 4 -). This energy saving stems primarily from the spectrum at low energies, 
below the magnetic excitation gap. This follows from the fact that the spin excitations are 
peaked around the wave vector Qaf at which /(q) is positive. Fig. |5] illustrates which 
part of the spectrum of Ihix(Qaf,w) increases/decreases AE X . This energy gain must be 
compared with the superconducting condensation energy AE C , which is the difference in 
internal energy between the (putative) normal and the superconducting state at T = O^-^l; 



Using the specific heat data shown in Fig. 6(a)^, we find the condensation energy to be 
?7 2.27 • 10 _4 meV per Ce. Compared to the high-T c cuprates where similar analyses have 
been performed^—, the considerably lower energy scales in the HF systems enable us to 




/(q) = I 1 [cos{q x a) + cos(q y a)} + ^{a, c, q) 



(1) 



AE c = U N (T = 0)-Us(T = 0). 



(2) 
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perform a quantitative analysis of the data in terms of an accessible putative normal state. 
As noted above, extrapolating the spin excitations from above T c is in good qualitative 
agreement with the excitations of the field-driven normal state at the lowest temperatures. 
Furthermore, the electronic specific heat of both the superconducting and normal state can 
be reliably determined since phonons do not contribute at such low temperatures. Despite 
this apparent advantage of HF systems, AE C and AE X have not received much attention 
in the context of HF SC. AE C has been determined for CeCoIn^, a compound whose 
proximity to quantum criticality is not yet certain, since SC sets in before AF order can 
develop. Our study represents the first determination of the savings in both the exchange 
energy and condensation energy for a superconductor near an AF QCP, as well as for any 
unconventional low-temperature superconductor. 

Our observation that the magnetic exchange energy saving is more than one order of 
magnitude larger than the condensation energy implies that AF excitations are the primary 
driving force for SC. A comparable factor of exchange energy saving over condensation en- 
ergy has recently been observed in the unconventional superconductor YbE^CusOe.^- As 
described above, the temperature dependence of a; ga p(T) in CeCu2Si2 follows a rescaled BCS 
form. For a conventional BCS superconductor, where ©d 3> w gap , the saving in potential 
energy is enhanced over the condensation energy by a factor that depends logarithmically 
on the ratio of Debye temperature ©d and superconducting gap w gap (T = 0)2£. The corre- 
sponding enhancement factor over the condensation energy in CeCu2Si2, where the Kondo 
temperature Xk ~ 15 K replaces the Debye temperature ©d, turns out to be two. The 
fact that the observed magnetic energy saving is more than a factor 20 larger than the 
condensation energy, indicates a large loss in kinetic energy. A natural origin for the latter 
lies in the Kondo effect, since the kinetic energy of the quasiparticles appears through the 
Kondo-interaction term. Because superconducting pairing in CeCu2Si2 occurs in the spin- 
singlet channel, the opening of the superconducting gap therefore weakens the Kondo-singlet 
formation and, by extension, reduces the spectral weight of the Kondo resonance^. 

Comparison with other unconventional superconductors 

Our understanding of the magnetic exchange energy saving in the HF superconductor 
CeCu2Si2 near its AF QCP naturally leads us to ask whether the effect is universal. SC- 
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induced enhancement of the spin-fluctuation spectrum in some frequency range has also been 
observed in the high-T c cuprates such as YBa2Cu307_<r£, iron pnictides such as K- or Co- 
doped BaFe 2 As2^>^ or FeTei-^Se^^, as well as two other HF compounds, UPd 2 Al^i2i 
and CeCoIn^. However, there are some striking differences between the spectrum observed 
in CeCu2Si2 and those seen in the other superconductors. In contrast to CeCu2Si2, where 
SC and long-range AF order exclude each other, SC in UPd2Ai3 occurs inside the anti- 
ferromagnetically ordered part of its magnetic phase diagram, which is far away from any 
QCP—. Whether a QCP underlies SC in the cuprates, the iron pnictides, or CeCoIn 5 , is yet 
to be established. The normal state magnetic response of S-type CeCu 2 Si 2 at Qaf slows 
down considerably, when lowering the temperature, indicating its proximity to a QCP, and 
displays pronounced dispersion. CeCu2Si2 represents, therefore, the only system in which 
we can unambiguously establish the linkage between AF quantum criticality and unconven- 
tional SC, even though the effect may well prove to be broadly relevant. In comparison to 
other HF superconductors, the inelastic spin response in CeCu2Si2 is broad in energy and 
extends beyond 10 times the gap value, while in CeCoIns a rather sharp, resolution-limited 
spin resonance is founcr^. Furthermore, unlike CeCoIns the temperature dependences of 
the spin excitation gap in CeCu2Si2 and UPd2Al3^ do follow the expected BCS form. In 
comparison to CeCu2Si2, UPd2Al3 also exhibits a dispersive spin excitation starting at the 
low-energy inelastic line (related to the edge of the spin gap>^) with a slightly higher in-plane 
mode velocity^. However, the situation in the cuprate superconductors is more complex, 
with an hour-glass like dispersion of the resonance mode^— . 

Experimentally the most prominent difference between CeCu2Si2 and other unconven- 
tional superconductors is the Q position of the spin excitation gap, which is observed in 
all reported unconventional superconductors at or close to simple commensurate positions 
with half- integer indices. E.g., in UPd2Al3 and CeCoIns it occurs at commensurate posi- 
tions, Q = (0 1/2) and (1/2 1/2 1/2) respectively^^. In contrast, in S-type CeCu 2 Si 2 
the gapped spin excitations are restricted to the vicinity of the ordering wave vector of the 
system, r m (0.215 0.215 0.53), which is incommensurate, far away from a simple commen- 
surate value. As a result, the opening of a spin gap becomes the major source of exchange 
energy saving. By extension, an additional excitonic resonance in x 5 (Q,o;) due to the su- 
perconducting stated would reduce the energy saving. This is a striking difference between 
CeCu2Si2 on the one hand, and CeCoIn^ and high-T c cuprate superconductors^ on the 
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other. 

In conclusion, our inelastic neutron scattering experiments in CeCu2Si2 reveal spin ex- 
citations associated with the AF (3D-SDW) QCP. These spin excitations are overdamped, 
dispersive and gapped in the superconducting state. Our quantitative estimate of both, the 
change in magnetic exchange energy and the superconducting condensation energy identifies 
the AF excitations as a major driving force for SC. AF QCPs are currently being explored 
in a variety of strongly correlated electron systems, including the new Fe pnictide super- 
conductors. Ba(Fei_ a; Co a ;)2As2r^, for instance, exhibits a T — x phase diagram very similar 
to the T — p phase diagram of CePd2Si2, raising the prospect that AF quantum critical 
excitations also drive the superconducting pairing in these new high-T c superconductors. 
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Methods 



High-resolution inelastic neutron scattering experiments were performed on the cold- 
neutron triple-axis spectrometer IN12 at the high-flux reactor of the Institut Laue-Langevin 
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in Grenoble/France. A vertical focusing graphite (002) monochromator and a doubly focused 
(vertical and horizontal) graphite (002) analyzer were used. The horizontal collimation was 
given by the neutron guide in front of the monochromator and 60' before the sample, while no 
collimation was inserted in the scattered beam. A liquid-nitrogen cooled Be filter was placed 
in the incident neutron beam to reduce higher-order contamination. The measurements were 
carried out with a fixed final wave vector kj = 1.15 A^ 1 which corresponds to a final neutron 
energy Ef = 2.74 meV and yields a high energy resolution AE « 57/xeV (FWHM, i.e., full 
width at half maximum). All experiments were performed on an S-type CeCu2Si2 single 
crystal (m « 2g). The crystal was mounted with the [110] axis vertical on a copper pin 
attached to the mixing chamber of a dilution refrigerator. The setup results in a [110] — [001] 
scattering plane. Data were taken at temperatures between T = 0.06 K and 1.7 K and 
in magnetic fields up to B = 2.5 T applied along the vertical [110] axis. The inelastic 
neutron scattering measurements were converted to units of /i|/(meV f.u.) by normalizing 
the intensities to the incoherent scattering of the sample. 



13 



Supplementary Material for "Magnetically driven superconductivity in 
CeCu 2 Si 2 " 

O. Stockert 1 , J. Arndt 1 , E. Faulhaber 2 3 , C. Geibel 1 , H. S. Jeevan 14 , 
S. Kirchner 1,6 , M. Loewenhaupt 2 , K. Schmalzl 5 , W. Schmidt 5 , Q. Si 7 , F. Steglich 1 

1 Max Planck Institute for Chemical Physics of Solids, Nothnitzer Str. 40, 01187 Dresden, Germany 

2 Institut fiir Festkorperphysik, Technische Universitat Dresden, 01062 Dresden, Germany 

3 Gemeinsame Forschergruppe Helmholtz-Zentrum Berlin - TU Dresden, 85747 Garching, Germany 

4 I. Physikalisches Institut, Universitat Gottingen, 37077 Gottingen, Germany 

5 Jiilich Centre for Neutron Science at Institut Laue-Langevin, 38042 Grenoble, France 

6 Max Planck Institute for the Physics of Complex Systems, Nothnitzer Str. 38, 01187 Dresden, 
Germany 

7 Department of Physics and Astronomy, Rice University, Houston, Texas 77005-1892, USA 
Sample characterization 

The present neutron scattering experiments were complemented by bulk measurements 
on the same S-type CeCu 2 Si 2 single crystal. The heat capacity was recorded using 
a compensated heat-pulse technique, while the ac susceptibility was measured with a 
homemade susceptibility setup and was recorded during the neutron scattering experiment 
simultaneously while the neutron data were taken. As shown in Fig. 6(a) the specific 
heat, plotted as C/T vs. T, exhibits a pronounced maximum indicating the onset of 
superconductivity at T c = 0.6 K. The shape of the anomaly at T c in comparison to heat 
capacity measurements on other CeCu2Si2 samples^ and its full suppression already in a 
magnetic field of B = 2 T clearly indicate that this is a transition into the superconducting 
state. The superconducting nature of this phase transition is further confirmed by ac 
susceptibility measurements displayed in Figs. Sl(b) and (c). A sharp drop of XaciT) 
at T c w 600 mK is observed with a large diamagnetic signal at lower temperatures. 
Furthermore, the susceptibility measurements yield an upper critical magnetic field B C 2 to 
kill superconductivity of B c2 — 1.7T for B || [lTO] at T = 0.07 K. 

From entropy the single-ion Kondo temperature Tk was deduced to almost coincide with 
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the lattice coherence temperature as obtained from the position of the low-temperature 
peak in the temperature dependence of the electrical resistivity^. Both are around 15 K. 
From the Doniach criterion together with the fact that our sample is almost quantum 
critical it follows that the magnetic energy scale Trkky ~ Tk ~ 15 K. 



Analysis of magnetic response 

The neutron intensity is directly proportional to the scattering function S(Q,oj). In 
the paramagnetic, normal state of S-type CeCu2Si2 the magnetic response was fitted by a 
quasielastic signal with Lorentzian lineshape (cf. Figs. 2(a) and (b) and Fig. 3(a)). Hence, 
the scattering function >S(Q, u) at momentum transfer hQ and energy transfer hu takes the 
form 

S qe ,mag(Q,u) = [n{Tvjj) + 1] • Imx(Q, u) = [n(flUj) + 1] 



vrr i + (frw/r) 2 ' 

at temperature T where n(ftw) + 1 = 1/[1 — exp(— hw/kT)] is the Bose factor and k is 
Boltzmann constant. Here, xo denotes the susceptibility and T the energy linewidth of 
the fluctuations (HWHM, half width at half maximum), which is inversely proportional to 
the lifetime r of the fluctuations. In the superconducting state the magnetic response was 
modeled using a modified scattering function taking the electronic density of states in a 
superconductor, Z(cu), into account, i.e., 



Z(u) 



= for u > u gap 



u2 ~~ Ugap 



otherwise. 



u g ap = ojg ap {T) is the value of the superconducting gap at a certain temperature T. 
Si n e,mag(Q, u>) = Sqe,mag(Q, • Z{oj) was then used to describe the data in the super- 
conducting state (cf. Figs. 2(a), (b) and 3(a)). In each case the scattering function was 
convolved with the instrumental resolution to fit the experimental data. 
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Determination of condensation and magnetic exchange energies 

The condensation energy AEc characterises the stability of the superconducting state 
(S) against a putative normal state (N) and is the difference in internal energy, or 

lim (G N (T, B = 0)- G S (T, B = 0)) = [ " M(T = 0) dB = \^VBl (S3) 

where Gs/Gn is the Gibbs free energy of the superconducting (S)/normal (N) state and B c 
is the thermodynamic critical field defined via \B 2 C = J^ c2 MdH. For CeCu2Si2, B c has 
been determined in Ref.-. However, the co-existence of small, magnetically ordered regions 
in the superconducting S-type CeCu2Si2 makes it necessary to perform all energy estimates 
on the same sample. Eq. flS3j) implies that AEc can be obtained from the specific heat data 
of Fig. 6(a): 

AE C = f C dT J T dT'^{C^{T')-C^T')) 

= r] 431 J/m 3 = rj 2.27 • 10™ 4 meV/Ce. (S4) 

Here and in what follows, we take the finite field (B = 2T) data as the putative normal 
state. From the lowest temperature (T < 0.1K) the data were extrapolated to T = and 
the contribution from the nuclear moments was subtracted. The contribution from below 
T = 0.1K to AE C is tiny (A££ <01K = 77 7.6 J/m 3 = r/4 ■ 10- 6 meV/Ce). The factor 77 > 1 
accounts for the fact that only the superconducting volume fraction (I/77) contributes to 
AEc, and it may very well be sample dependent. 

We now turn to estimating the difference in exchange energy between the normal and 
superconducting states, in order to ascertain whether the magnetic excitations contribute 
significantly to the condensation energy. Heavy-fermion metals are best described by the 
Anderson lattice model^. To estimate the magnetic exchange energy in such a system, it 
is sufficient to consider the magnetic limit of the Anderson lattice model, i.e., the Kondo 
lattice model: 

HKL = Yl C k C L C ka + J K S i ■ S c( r *)> ( S5 ) 
k,cr i 

where Sj is the localised moment at a cerium site which is coupled via the Kondo coupling 
Jk to the spd conduction electron spin density s c (rj) at the cerium site, and describes 
the conduction electron bandstructure. The Kondo model describes not only the formation 



16 



of heavy quasiparticles but also the AF phase and the quantum critical point, since Hkl 
implicitly contains the RKKY interaction among the Ce moments: 

Hrkky = / • Sj, (S6) 



i<j 



where the exchange constants Ijj ~ J^. 

We model the exchange interaction between the localised Ce-moments by including near- 
est neighbour and next nearest neighbour terms appropriate for the tetragonal, body-centred 
unit cell: 

/(q) = h[cos(q x a) + cos(q y a)] + I 2 [cos(q x a/2 + q y a/2 + q z c/2) 

+ cos(-q x a/2 + q y a/2 + q z c/2) + cos(-g x .a/2 - q y a/2 + q z c/2) 

+ cos(g x a/2 — q y a/2 + q z c/2) + cos(p x a + p y a) + cos(p y a — p x a)~\ , (S7) 

where a = 4.1 A and c = 9.9 A are the lattice constants for the a- and c-axis, respectively 
and ii and I2 are the nearest and next nearest neighbour exchange interactions. The energy 
saving in magnetic exchange energy of the superconducting ground state compared to the 
putative normal ground state is then given by 

AE X = E» -E* = ^J™ d -^± [n(M + 1] (S8) 

/tv I a pir/a ptv/c 

dq x / dq y J dq z I(q x ,q y ,q z ) 
-TV /a J —ii /a J — 7r/c 

I pir/a rTifa pn/c 

x Im X N (lx,q y ,qz,f^) - X S (<lx,qy,qz,?Uj) II dq x dq y dq z . 

I J ~ir/a J —Tr/a J —n/c 

Crystalline-electric-field effects split the J = 5/2 states of the Ce 3+ ion up into a ground-state 
doublet and a quasi-quartet at high energies (> 30 meV) and result in g-f actors g z ~ g± ~ 2, 
and an almost isotropic spin susceptibility^. A is a constant given byA = ?7-8-3/2 resulting 
from the eight symmetry equivalent incommensurate AF wave vectors in the first Brillouin 
zone and the fact that neutrons only detect moments and spin fluctuations perpendicular 
to the actual momentum transfer. Our sign convention in Eq. (1S3I) implies that a positive 
AE X is equivalent to an energy saving in the superconducting state. Different energy ranges 
will in general contribute differently in either decreasing or increasing the exchange energy 
as the system goes from the (putative) T=0 normal state to the superconducting state. 
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The exchange constants I\ and I2 of Eq. flST|) follow from three independent relations: 
(i) the observed dispersion of Fig. 7(a), (ii) the fact that S-type CeCu2Si2 is close to 
quantum criticality and meets the Doniach criterion and (iii) the functional form of the 
RKKY interaction. To estimate the magnitude of the exchange interaction Ii, we chose 
v exc from the observed dispersion of the overdamped excitations (see inset of Fig. 4(a)). 
Using the mean field expression of v exc for a three-dimensional, cubic lattice with I2 = 
0, we find I\ = 0.63 meV. Note, that an I2 > will lead to an increase in I\ for fixed 
dispersion, making I\ = 0.63 meV an estimate from below. Alternatively, estimating I\ 
via the Doniach criterion, we find that 7i « 0.6 meV when using Tk ~ 15 K and the fact 
that S-type CeCu2Si2 is almost quantum critical. The ratio I1/I2 is estimated from the 
distance-dependence of the RKKY interaction. This interaction shows oscillatory behaviour 
with an envelope that falls off as a function of the inverse of the distance r between the 
Ce moments. For free electrons, this function is a simple power law and the period of 
the oscillating function is set by 2kpr. In CeCu2Si2, the RKKY interaction is mediated 
by the spd conduction electrons. For a reliable estimate of the ratio in CeCu2Si2 

we combine band structure calculations for the non- magnetic La-homologue LaCu2Si2 of 
reference^ with the RKKY interaction for non-spherical Fermi surfaces obtained by L. Roth 
et al.-. As found by Roth et al., the RKKY interaction shows a 1/r 3 dependence even for 
non-spherical Fermi surfaces except for directions where the effective band mass diverges. 
We therefore set I1/I2 ~ t\ cos(cj)i) / (rf cos(02)) or J 2 ~ 0.35ii cos(02)/ cos(0i), where r\ 
ir2) is the distance between nearest (next nearest) neighbour Ce atoms and cos(0i) and 
cos(02) are oscillating factors which not only depend on the magnitude of rxjr^ but also 
on the direction for non-spherical Fermi surfaces. The oscillating factors are determined by 
the wave vector differences of calipering pairs of points of the Fermi surface along the [110] 
for the next nearest neighbour in the basal plane and [111] directions for the next nearest 
neighbour in the unit cell centre^. 

As a result, the next nearest neighbour exchange constant /Jp asa l = — O.O8/1 for two Ce 
moments in the basal plane differs from the exchange constant /tenter _ q gQ^ f or nex t 
nearest neighbours in the unit cell centres. We are only interested in a lower estimate for 
the exchange energy gain and use a simple average I2 = (j| ,asa ^ + 2/2 entre )/3 ~ 0.38/i for 
the next nearest neighbour exchange constant in Eq. (IS7j) . Note that there are twice as 
many next nearest neighbours in the unit cell centre than in the basal plane. 



18 



The volume fraction 77 enters Eq. (1S8I) . since the magnetically ordered regions in S and N 
yield (essentially) identical responses for B = and B = 2T. Therefore, the actual change 
in exchange energy between S and N is larger by a factor rj. As a result, the ratio between 
AEc and AE X will be independent of 77. Nonetheless, an estimate of the volume fraction 77 of 
our sample can be obtained from the ordered moment associated with the elastic magnetic 
response and the weak anomaly in the heat capacity yielding m 0.02 /ie/Ce. Comparing this 
to the ordered moment in the A-phase, fi or d ~ 0.1 /is/Ce yields rj ~ 1.25. 

The spin susceptibility in the normal state has been parametrized as 

W(Q,«) - 1 + e(Q % AF) ^/b (S9) 

Xo 2tu/T 

1 + {£/a)\a{Q x - Ql ¥ )f + {i/af{a{Q y - Q V AF )) 2 + (£/cf(c(Q z - Q AF )) 2 1 + 4(^/T) 2 ' 

In the superconducting state, the susceptibility has been modeled as 

OJ 

^gap 

for oj > ujgap and zero otherwise. These expressions model very well the experimental data 



Im X 5 (Q,a;) = w ^(Q,^), (S10) 

w — CJ„ 



for Imx^Q,^;) and Imx s (Q,a;) in the vicinity of Qaf at small energy transfers where u g 
is independent of Q. Away from Qaf and at energy transfers above 0.5 meV Imx^Q,^;) 
and Imx 5 (Q,u;) are identical (within the experimental error bars) and therefore do not 
contribute to AE X . 

Putting everything together, we obtain AE X = ?74.8-10 _3 meV/Ce and AE X / AEq = 21.1. 
The parameters in Eq. (IS9l) and (ISlOj) were obtained from fits to the experimental data at 
T = 0.07 K in the superconducting (B = 0) and normal (B = 2T) state as £ = 20 — 25 A, 
T N = 0.11 meV, Xo = 15.64/4, r s = 0.225 meV, and Xo = 8.69/4- We checked that the 
associated static structure factors integrated over the full Brillouin zone ( "local moment sum 
rule") in the normal and superconducting state yield identical results (within a 5% error). 
The (in reciprocal space) isotropic fit to Imx^ (Imx s ), see^, suggests that the value for AE X 
is an estimate from below. We also checked that even if J 2 were zero, AE X would still be 
one order of magnitude larger than AEc- 

A realistic modelling of the dynamic susceptibility of the normal and superconducting 
states in the entire momentum and frequency range should include the overdamped, disper- 
sive excitations. This can be accomplished by parameterising the normal state susceptibility 
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by the SDW form 

which includes an overdamped excitation at 

u = u exc = ±b~ V^" 2 + (Q - Qaf) 2 (S12) 

where £ is the correlation length and 6 = u^, follows from the dispersion relation of Fig. 
4(a). Eq. ( IS 11 1) reproduces the observed Lorentzian lineshape of the overdamped, dispersive 
mode for a fixed u and Q ~ Q exc (where Q exc is a solution to Eq. (IS12P ): 

lmx N (Q,u) = ^ZE _ (S13) 

7T(Q-Q exc )2 +u; 2 /r2 ' 

where T = 2T\Qaf — Qexc I and x = 7t/(2£ 2 |Qa,f — QexcD- Fitting the Q-dependence 
of the data at different energies gives rise to the same width, as required by Eq. (IS 1 3[) . 



Imx^Q,^) and Imx 5 (Q,oj) are identical at energy transfers above 0.5 meV and are still 
sizeable in magnitude for Q sufficiently away from Qaf- We found that the best fit to the 
superconducting lmx s could be obtained from 



W(Q^) = &(io gap -Lo) 3 _ U 3 lm X N (Q^). (S14) 

^ gap) 

Taking Xo = 7.0 T N = 25meVA 2 and %o = 6.9 V s = 25meVA 2 reproduces all 
features of our data (cf. Fig. [7]) and results in a saving in exchange energy AE X = t] 4.7 ■ 
10~ 3 meV/Ce. The slightly modified parameters X q = 10.0 /ig, T N = 20.0 me VA 2 and 
Xo — 10.0 T s = 20.0 me VA 2 capture the broad features in Imx(Q,w) and result in a 
fit of somewhat lesser quality to the data, but still yield AE X = t] 4.1 • 10 _3 meV/Ce. We 
therefore conclude that AE X ^> AEc is a stable observation insensitive to the details of the 
fitting and results primarily from the changes of Im%(Q,o;) at low energy transfers in the 
vicinity of Qaf as superconductivity sets in. The saving in exchange energy is more than 
an order of magnitude larger than the condensation energy, thus identifying the build up 
of magnetic correlations near the AF QCP as the major driving force for SC in CeCu2Si2- 
It is important to note, that an increase in AE X comes from the opening of the spin gap 
and not the 'resonance '-like feature above the spin gap, which tends to reduce the energy 
saving as illustrated in Fig. 5. This figure shows the difference of Im X s and Im x ■ The 
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blue area marked with a '+' contributes to an increase in AE X whereas the green area 
leads to a reduction in AE X , as a consequence of /(Qaf) > and similar for the other 
wave vectors as follows from Eq. (lS9p . For e.g. CeCoIn^&ii, the shifted spectral weight 
in the superconducting state has to be such that the resulting positive (since /(Q) > 0) 
green area exceeds the blue area for an overall saving in exchange energy. Therefore, the 
spin resonance at low energies observed in CeCoIns and in the cuprates contributes to an 
increase in AE X . A sharp resonance occurs in response to the superconducting state in 
predominantly two-dimensional superconductors and adds to a saving in exchange energy. 
It is however not expected in CeCu2Si2, which is a 3D superconductor as deduced from the 
nearly isotropic upper critical field^. Instead, as shown in e.g. Figure 2(a), the magnetic 
response in CeCu2Si2 is broad and extends to more than ten times the gap energy in contrast 
to e.g. CeCoIns or the cuprates where a sharp spin resonance has been observed. 

Our results imply that there is a sizeable "kinetic" energy loss in CeCu2Si2- As described 
in the main text, superconductivity in CeCu2Si2 occurs in the spin-singlet channel. As a 
result of the opening of the superconducting gap, the Kondo-singlet formation is weakened 
and the spectral weight of the Kondo resonance is reduced. The spectral weight sum rule 
dictates that the integration of the single-electron density of states over all energies is un- 
changed (and is equal to one). As a result, the lost spectral weight from low energies (within 
the scale of the Kondo temperature) must be transferred to higher energies, on the order 
of the on-site Coulomb interactions among the /-electrons where the incoherent /-electron 
excitations reside. This energy loss in a heavy fermion superconductor should therefore be 
distinguished from the ordinary kinetic energy loss of a classical superconductor—. 
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Figures 




FIG. 1: Schematic phase diagram around the QCP, crystal structure and nearly 
isotropic spin fluctuations of CeCu2Si2. (a) Schematic T — g phase diagram of CeCu2Si2 
in the vicinity of the quantum critical point (QCP) where the antiferromagnetic (AF) phase van- 
ishes as function of the effective coupling constant g. Superconductivity (SC) is observed around 
the QCP and extends far into the paramagnetic (PM) regime. Composition as well as hydro- 
static pressure can be used to change the coupling constant g and to tune the system to the QCP. 
The positions of the A-type and the S-type single crystals in the phase diagram are marked, (b) 
Tetragonal crystal structure (space group: IA/mmm) of CeCu2Si2- The nearest and next-nearest 
neighbour interactions between the cerium atoms are labelled by I\ and I2. It should be noted 
that the distances between next-nearest neighbour Ce atoms in the basal plane and out-of-plane 
are almost identical, (c) The spin fluctuations at T = 0.06 K and B = and at an energy transfer 
tvjj = 0.2 meV. The anisotropy factor between the [110] and the [001] directions is about 1.5. Note 
that the correct aspect ratio [110]* : [001]* has been taken into account although the axes are 
labelled in reciprocal lattice units (rlu). Black dots mark the (Q,w) positions data were taken. 
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FIG. 2: Inelastic neutron scattering spectra in the normal and superconducting states 
of CeCu2Si2- Energy scans (neutron intensity S = S e i a + S qe / inejmag versus energy transfer Hlj) in 
S-type CeCu2Si2 at Q = Qaf = (0.215 0.215 1.458) in (a) the superconducting state at T = 0.07 K, 
B = and in (b) the normal state at T = 0.8 and 1.7 K, B = and T = 0.07 K, B = 2T. For 
comparison the magnetic response at an arbitrary, general Q position Q = Q ar b = (0.1 0.1 1.6) at 
T = 0.07 K, B = is also plotted in (a). The inset in (a) shows the magnetic response at Qaf 
(T = 0.07 K, B = 0) extending beyond fko = 2meV. Solid lines represent fits to the data comprised 
of the incoherent and coherent elastic signal S e i a and the quasielastic/gapped inelastic magnetic 
response S qe /i ne ^ rnag (dashed lines) with Lorentzian lineshape convolved with the resolution. The 
gapped magnetic response at T = 0.07 K, B = has been modelled by a quasielastic Lorentzian 
line taking into account a spin gap with a value hu} gap and an enhanced density of states above 
the gap as for the electronic gap of a BCS superconductor. The error bars represent the statistical 
error. 
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FIG. 3: Magnetic response, relaxation rate and spin gap at the AF wave vector of 
CeCu2Si2. (a) Magnetic response S qe /i ne ^ mag at Qaf and T = 0.07 K in the superconducting 
and the normal state, i.e., at B = and B = 2T, as extracted from the data displayed in Fig. 2. 
The data have been put on an absolute intensity scaled. Below huj ~ 0.1 meV (~ 1.6x the 
instrumental resolution) the errors in S ma g increase strongly (some data points fall even outside 
the plotted range), since the strong elastic scattering S e [ a is subtracted from the total scattering 
to receive S mag and due to small uncertainties in the resolution function. These uncertainties are 
the same for both data sets and do not play a role since only the difference is analysed for the 
estimation of the exchange energy saving, (b) Linewidth V vs. temperature T of the quasielastic 
magnetic response at Qaf in the normal state as yielded by fits to the data shown in Fig. 2. 
Plotted here is T(T) - T Q vs. T, with T = 0.112meV. The solid line F(T) -T = aT 3 / 2 (with 
a = 0.061 meV/K 15 ) is the expected behaviour near a 3D SDW QCP. (c) Temperature dependence 
of the spin excitation gap foui gap at Qaf together with the scaled d-wave BCS superconducting gap 
function (solid line). The error bars denote the statistical error. 



26 




FIG. 4: Dispersion of the magnetic response in superconducting CeCu2Si2. (a) Wave 
vector Q dependence of the magnetic response around Qaf in S-type CeCu2Si2 in the supercon- 
ducting state at T = 0.06 K for different energy transfers frw. The scans are shifted by 100 counts 
with respect to each other. Solid lines denote fits of two peaks with Gaussian lineshape to the 
data, while dashed lines are only guides to the eye. From the linewidth at small energy transfers 
a dynamic correlation length £ ~ 25 A is inferred. Inset: Dispersion of the magnetic excitation 
around Qaf at T = 0.06 K as a result of the fits to the Q scans. The solid line indicates a fit 
to the data with a linear dispersion relation yielding a velocity v exc = (4.44 ± 0.86) meVA (for 
comparison spin-wave velocities in other HF metals, UPd2Al3: v = 10 — 15meVA2£, URu2Si2: 
v ~ 45 me V A— ). The error bars represent the statistical error, (b) Colour coded intensity plot of 
the data displayed in (a) and additional data, clearly indicating the dispersion of the gapped spin 
excitation. Black dots mark the (Q,w) positions data were taken. 
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FIG. 5: Schematic plot of the imaginary part of the dynamic spin susceptibility 
Imx(QAF>^) hi the normal (N) and superconducting (S) states. The dynamic correla- 
tion function 5(Qaf 5 w ) of Figure 3(a) is related to Imx(QAFi^) via the fluctuation-dissipation 
theorem after de-convolving the data with the instrument's energy resolution function. The blue 
area marked with a '+' contributes to an increase in AE X whereas the green area (marked with a 
'-') leads to a decrease in AE X . The fact that the opening of the gap contributes to the saving in 
exchange energy is a consequence of /(Qaf) > at the wave vector Qaf, where Imx N ' S (Q, ui) is 
peaked. 
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FIG. 6: (a) Specific heat C of the same S-type CeCu2Si2 single crystal which was studied by 
neutron scattering, plotted as C/T versus temperature T in zero magnetic field B = and a 
magnetic field B = 2T applied along [110]. (b) Temperature dependence of the ac susceptibility 
Xac at B = 0. (c) Magnetic field dependence of Xac at a temperature T = 0.07 K and in a magnetic 
field B || [110]. The data with large error bars were taken when sweeping the magnetic field while 
all other data were measured during neutron scattering scans where the magnetic field was constant 
for a long time. 
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FIG. 7: Magnetic response at Qaf i n (a) the normal and (b) the superconducting states of CeCu2Si2 
(same data as in Fig. 3(a)). Solid lines are fits to the data (see text). For energies Two > 0.5 meV, 
the dynamical susceptibility of the normal and superconducting states coincide (within the exper- 
imental error bars) and thus do not contribute to the exchange energy saving. 
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